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Plasma instabilities excited in collisionless shocks are responsible for particle acceleration, gen- 
eration of magnetic fields, and associated radiation. We have investigated the particle acceleration 
and shock structure associated with an unmagnetized relativistic jet propagating into an unmag- 
netized plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are 
swept up to create a partially developed hydrodynamic-like shock structure. The shock structure 
depends on the composition of the jet and ambient plasma (electron-positron or electron-ions). 
Strong electromagnetic fields are generated in the reverse, jet shock and provide an emission 
site. These magnetic fields contribute to the electron’s transverse deflection behind the shock. 
We have calculated, self-consistently, the radiation from electrons accelerated in the turbulent 
magnetic fields. We found that the synthetic spectra depend on the Lorentz factor of the jet, its 
thermal temperature and strength of the generated magnetic fields. The detailed properties of the 
radiation are important for understanding the complex time evolution and/or spectral structure in 
gamma-ray bursts, relativistic jet shocks, and supernova remnants. 
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1. Introduction and Simulations 


Particle-in-cell (PIC) simulations can reveal the physical mechanism of particle acceleration 
that occurs in the complicated dynamics within relativistic shocks. Recent PIC simulations of 
relativistic electron-ion and electron-positron jets injected into an ambient plasma show that accel- 
eration occurs within the downstream jet [4]. 

The effects of different mass ratios (mi/m e = 1,20) are shown in Figure 1 [2]. Figure l(b, and 
d) are from an electron-positron simulation as compared to an ion-electron simulation shown in 
Figure l(a, and c). Forward and reverse shock structure is generated in both cases but differences 
in the shape of the generated shock structure and in the magnetic fields are immediately apparent. 
In the electron-ion case a sharp contact discontinuity (CD) with strongly enhanced magnetic field is 
produced (Figs, la and lc). In the electron-positron case the maximum associated with the forward 
shock is significantly increased. In this case, both positrons and electrons in the ambient appear 
swept up by the front of the jet without the strong interactions seen in the electron-ion case, and the 
jet density exhibits a small peak associated with the shocked ambient particles. For both electron- 
ion and pair plasmas cases magnetic fields shown in Fig 1 (c, and d) increase in the reverse shock 
region, but the maximum magnetic field occurs at the contact discontinuity (CD) in the electron-ion 
case. Magnetic (Figs, lc and Id) and electric field (not shown) in the transverse directions strongly 
fluctuate in this region. 
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Figure 1: Snap shots of ion density (a), positron density (b), and magnetic field strength. If ((c) for electron- 
ion plasma, and (d) for pair plasma) at = 375. 
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2. Synthetic spectra from simulations 

We have calculated the radiation spectra directly from our simulations by integrating the ex- 
pression for the retarded power, derived from the Lienard-Wiechert potentials for a large number 
of representative particles in the PIC representation of the plasma[5]). In order to obtain the spec- 
trum of the synchrotron/jitter emission, we consider an ensemble of electrons selected in the region 
where the Weibel instability has fully grown and where the electrons are accelerated in the self- 
consistently generated magnetic fields. 

Figure 2 shows how our synthetic spectrum matches with spectra obtained from Fermi obser- 
vations. Figure 2a shows the observed spectra in vF v as modeled by[l] at five different times. The 
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Figure 2: Comparison of synthetic spectra (b) with Fermi observations for GRB 080916C (a). Figure 2a 
shows modeled Fermi spectra in vF v units at five different times. A flat spectrum would indicate equal 
energy per decade in photon energy. The changing shapes show the evolution of the spectrum over time. 
Figure 2b shows PIC generated spectra for jet Lorentz factors of y = 10, 20, 50,100, and 300. For each case 
we show two simulations represented by the same color for cold (thin, lower lines) and warm (thick, ipper 
lines) electrons. The low frequency slope is approximately 1. 

red lines in Fig. 2a and 2b indicate a slope of one, and except for the Fermi spectrum at time “a” the 
low frequency slopes are all approximately one. As shown in Fig. 2a, prompt emission spectra of 
GRB 0809 16C exhibit a remarkably constant shape (low-energy and hi-energy slopes are different 
but do not evolve significantly). On the other hand the peaks of the SED vary significantly with 
time. Based on our PIC simulations, this behavior can be interpreted in terms of declining bulk 
Lorentz factor. 

Behind the reverse shock electrons are accelerated and strong magnetic fields are generated as 
shown in Fig. [4, 2]. Therefore, this region would seem more likely to produce the emission that is 
observed. In the future we will examine PIC computed spectral changes over time varying plasma 
conditions such as plasma composition and ambient magnetic fields along with jet Lorentz factors 
and jet thermal temperatures. We are currently simulating much larger systems to obtain synthetic 
spectra at different times in the nonlinear development where different particle acceleration rates 
and magnetic field strengths apply. 


3 


Radiation from shock-accelerated particles 


3. Discussion 

A double shock structure (bow and jet shocks separated by a contact discontinuity region) is 
formed and electrons can be accelerated due to the Weibel instability in both shocks. Since we 
calculate the radiation from the electrons in the observer frame calculated spectra can be compared 
directly with observations. As shown in Fig. 1, the strongest electron acceleration and strongest 
magnetic fields are generated in the reverse, jet shock. Therefore, in this simulation this region 
would produce the emission that is observed. 

Emission obtained using the method described above is obtained self-consistently, and auto- 
matically accounts for magnetic field structures on the small scales responsible for jitter emission. 
By performing such calculations for simulations using different parameters, we can investigate 
and compare the different regimes of jitter- and synchrotron-type emission[3]. This addresses the 
observed low frequency GRB spectral index violation of the synchrotron spectrum line of death[3]. 

More extended simulations are required for detailed comparisons with GRB spectra. In par- 
ticular - , the nonlinear development of the Weibel instability must be followed for a sufficiently 
long time, and radiation spectra must be sampled for the various plasma conditions suggested by 
dynamic models of jet evolution. 
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